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'Abstract 

Events of giant conductance and anomalies of the phase transmission time for holes, are theoretically investigated 
within the multicomponent scattering approach. Based on this model, new analytical expressions for unitarity re- 
lations in the uncoupled hole transport are obtained and directly applied to study the behavior of the conductance 
and the phase transmission time in a double barrier resonant tunneling (DBRT) and a superlattice GaAs-cladding 
layer/(74M.s/GaA.s)"/GaAs-cladding layer. Clear-signature evidences of giant conductance phenomena for hole trans- 
mission without valence-band mixing through a DBRT and a superlattice were found. The giant conductance effect 
losses robustness by manipulating the number of superlattice layers and by including the valence-band particles coupling 
as well. Phase time through the heterostructure exhibits extremal dependencies in the gaps and in the barriers, as those 
reported before for electrons. We have detected an earlier arrival phase time for the propagation of both flavors of holes 
within the barrier, in the order of few tenths of picoseconds. An appealing filter-like effect is presented, whenever a 
Selective confinement strength arises independently for both flavors of holes in the uncoupled regime. Our results also 
.prescribe noticeable evidences for both uncoupled and coupled hole fluxes, similar to those foretold by Hartman, upon 
transmission of electrons through opaque barriers. 

J^eywords: tunneling, phase time, elastic quantum scattering 
■PACS: 73.43.Jn, 34.80.Bm 



1. Introduction 

The impressive development of low-dimensional elec- 
tronic and optoelectronic devices, brought a new urgency 
.to the essential measurement and modeling of charge car- 
'riers transmission time and particularly of the tunneling 
time through specific potential regions. When electrons 
and holes are involved, the low-dimensional device response 
idepends on the slower charge-carrier's traveling time 
'Undoubtedly, for such systems it is crucial to study the 
.tunneling time of holes. An equally important reason for 
'address investigations on heavy holes (hh) and light holes 
,{lh), is their intrinsic band mixing effects. Since the hole's 
transport is essentially an interconnected-multichannel pro- 
cess, the theoretical calculation of transmission quantities 
is in some sense more cumbersome than that for electrons. 

Currently, the research of under-synchronic anomalous 
'scattering events, is one of the major areas in the quantum 
transport theory 0-01 • Regarding the long-standing con- 
troversy around superluminal velocities and total trans- 
mission inconsistency^together with their wide- known ano- 
malous character [§| , they are welcome due to all pos- 
sible benefits they open up for design and development of 
nano-optoelectronic devices [l2. 13 1. There is a large ac- 
cumulated knowledge about these topics -mostly dedicated 
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to electrons, optical pulses and electromagnetic waves 
lol- 15|-. accurately enough to unambiguously provide solid 
grounds for superluminal propagation, the paradoxical Hart- 
man premonition and legitimacy of the phase time ap- 
proach. Worthwhile noticing however that, hole transport 
analogous phenomena has been less studied yet. 

Resonant mechanisms between propagating modes os- 
cillations and the self oscillations of the system are the 
fundamental cause, which gives rise to an isolated maxi- 
mums of transmission, often referred to as giant conduc- 
tance (GG) phenomena. Events of GG have been observed 
for electrons and optical pulses in many finite periodic sys- 
tems, which possess time reversal symmetry, namely: bal- 
listic mesoscopic conductors, luminal pulses and super lat- 
tices [i3,[ii-[22|. Recent theoretical studies in molecular 
devices suggested that could rise up particularly large val- 
ues for dynamic conductance j23| and novel experimental 
achievements pointed out to a giant enhancement of elec- 
tronic tunneling in a two-dimensional lattice of graphenejl^, 
[25! . Lately, was found a strong dependence of the conduc- 
tance fluctuations with the density of states in graphene 
nanoribbons, which represent an alternative i nsig ht to study 
the band structure of these physical systems [26|. 

Owing to its relevance for quantum transport, the tun- 
neling time has played a key role for a correct evaluation 
of electronic devices A large number of proposals in 
reply the challenge to settle how long does it takes to a 
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particle to tunnel a single potential barrier [30| , have been 
put forward d, 0, [s^l . The lively and long standing de- 
bate concerning this matter d M, [H, II^^Hi -presumed 
not to finish completely-, it seems to have finally found 
its cornerstone in precise experimental measurements for 
photons and optical pulses [l2l . 13 . 36 -3^ 
uine measurements put to the 



These gen- 
all else, the 



test, before 

very question loosely formulated by McCoU and directly 
bear on truly temporal features of tunneling phenomenon. 
Besides, they are consistent with the phase time concep- 
tion derived from group delay within the stationary-phase 
method (28l . [29I . An additional motivation was the 

accurate description of above mentioned experiments by 
careful phase-time calculations [l^, U | . 

In this paper we will present numerical calculations of 
the conductance and the phase time for an arbitrary num- 
ber of alternate layers of /// — V materials, in the frame- 
work of the multicomponent scattering approach (MSA)[^. 
The main idea is to search evidences of anomalous events 
in hh and Ih transmission. In Secj^jsome theoretical re- 
marks within the MSA are exposed. Brief complemen- 
tary insights to the uncoupled regime formula of scatter- 
ing matrix unitarity relations and phase time will also be 
presented. Analytical expressions for the phase time and 
the conductance are applied in Sec|3]to study the anoma- 
lous events associated with these physical quantities in the 
double barrier resonant tunneling (DBRT) and the super- 
lattice. Phase time behavior for holes shows a very good 
agreement with the previous results for electron-tunneling 
through superlattices[ll|. Evidences of the giant conduc- 
tance phenomena for uncoupled hole transport are shown. 
Finally, in Sec|3]we summarize and conclude. 

2. Conduction and phase time in the multicompo- 
nent scattering approach 

The MSA is based on standard (N x N) k • p effective 
Hamiltonians and the multichannel transfer matrix (TM) 
method. In this approach the scattering amplitudes can be 
obtained straightforwardly The main purpose here, is 
to evaluate quantities of synchronic multi-mode transport 
for massive charge carriers moving through multi-layered 
structures at coupled and decoupled channel regimes. 

Owing to consistency, we found useful to recall few 
basic outlines of the Kohn-Liittinger (KL) two-band model 
(39I . due to its widely accepted accuracy for describing 
dynamics of elementary excitations as well as electronic 
properties in the valence band, within the framework of the 
multiband effective mass approach. This model, is highly 
common regarding to the field of hole quantum transport 
through periodic heterostructures. It considers in the k p 
approximation the highest two valence bands degenerated 
in the F-point of the Brillouin Zone. The usual (4 x 4) 
KL Hamiltonian [39] in the Broido-Sham representation 
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Figure 1: Schematic illustration of multichannel hh 

and Ih quantum transport processes through a GaAs — 

c\a.dding \aycT:(V)/{Al As / GaAs)" /GaAs — cladding layer{R) 
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being 



Ai.2 = (71 ±72) 
Bi,2 = (71 ± 272), 



in atomic units. The matrix elements Hij show a direct 
dependency of the semi-empirical parameters of Liittinger 
7i (i = 1,2,3), that typify each slab of the heterostruc- 
ture. While V{z) and Kt, stand for the periodic poten- 
tial, and the two-components in-plane gMasi-momentum, 
respectively Q. 

To study the quantum transport properties of holes 
through semiconductor heterostructure, we consider the 
system shown in Fig.([T]), described by the Kohn-Liittinger 
Hamiltonian in the expression ([1}. The four accessible 
channels of the system, denoted as: channel 1 : hh 
channel 2 : lh_ 



1/2: 



;m, aenotea as: cnannei i : /i/i+3/2, 
channel 3 : lh^i/2, and channel 4 : 
hh_3/2i are considered simultaneously, following one of the 
bare bones of the MSA. We emphasize this upmost advan- 
tage, as it allows clearly identify inter- and intra-subband 
hole transitions together with transitions from Kramer- 
up toward Kramer-down states to time reversal, and the 
other way around. The latter is not the usual case. In 
this approach, the segment between z^ and Z3 in Fig. ([T]), 
is defined as a cell and represents the period of the het- 
erostructure superlattice (3|. 

Combining the (4 x 4) Kohn-Liittinger Hamiltonian 
with the TM formalism, a direct connection with the scat- 
tering theory is established Q, then we have Afs„(zfl, zl) = 



2 



Af~^ M fd{zR, zl)./^- Here. J\f defines a relation between 
the traveling wave vectors of the tunneling problem and 
their derivatives, with the A^— dimensional state vectors 
that describe a propagating or evanescent mode. This 
crucial transformation of the MSA, connects two differ- 
ent kinds of TM. The first class of TM, M/d(z_R, zl), re- 
lates the linearly independent state functions and their 
derivatives {fd) between two points of the heterostructure 
in the form M fd{zu, zl) = N{zji) ■ N{zl)~^, where the 
matrix N{zi) can be constructed once the linearly inde- 
pendent solutions are obtained The second kind of 
TM, M sv{zr, zl), whose elements describe a propagat- 
ing or evanescent mode depending on the energy, connects 
the state vectors {sv) in two different points of the sys- 
tem. The mathematical structure of these two matrices 
is conditioned by requirements of time-reversal invariance 
and flux conservation f4l| . A very strong criterium to ver- 
ify the calculations is based in the determinant of these 
transfer matrices which must be constant and equal to the 
unity, mostly. 

Then applying the well known relations between the 
a. (3 
7 S 
r t' \ I — I 

^ ^, j blocks [41|, one can directly evaluate the 
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and the scattering matrix 



relevant physical quantities of the scattered transport the- 
ory. Initially we have to obtain the complex transmission 
amplitudes 3 

t = a. f3S^^^, 



(2) 



for each transition from the j-th incoming channel (in- 
cident propagating mode) to the i-th outgoing channel 
(transmitted or reflected propagating modes). It is worth 
noticing that within the MSA, none of the incoming ampli- 
tudes is required to cancel, as it is commonly done. Hence, 
an authentic multichannel and synchronic description of 
the hole transmission processes is naturally achieved. This 
is to say that, we are able to distinguish accurately the 
contribution from each input channel to the scattering pro- 
cess, modulated by the simultaneous presence of the rest 
of the channels, which are, in principle, equally accessible 
to the incident hole stream. None of the usual simplifica- 
tions of the original fc • p Hamiltonians are invoked, thus 
the time reversal symmetry is assured. Given the trans- 
mission amplitudes t^ , it is possible to evaluate the total 
transmission probability to channel i 



G, = iey7Th)J2\Uj\\ 



(3) 



the total two-probe Landauer conductance defined as 

G=(e2M)Efel'' (4) 



and the transmission amplitude phase for each of the 16 
available paths ij Q, 9ij = arctan{ti../tfj..}. The trans- 



mission phase time is obtained from Q 



(5) 



whose advantages and applicability criteria has been widely 
discussed in the literature [3, [l^, U , l3) 13 • The sums in 
equations ^ and (|1]) run over all the open channels and 
the expressions ([31)- ([S]) and ([S]) are defined for the n cells 
periodic structure showed in Fig.([T|). Hereafter E stands 
for the incoming particle's energy. 

Tunable parameter Kt, is essential to modify the in- 
plane dynamics of holes, i.e. the degree of freedom trans- 
verse to the main transport direction, also refereed as the 
valence-band mixing for holes. This later is crucially rele- 
vant in order to understand the main problem envisioned 
here, and how we face it within the MSA. Whenever the 
gMfflsi-momentum parallel to the interfaces, Kt ~ 0, the 
crossed transitions are not allowed and the system is said 
to yield an uncoupled hole regime. Under this particu- 
lar condition, only direct incoming-outgoing channels are 
connected, therefore (4 x 4) blocks a, (3, 7 and S of the 
TM Msv are diagonal matrices. The later, straightfor- 
wardly leads us to new unitarity relationships, that satisfy 
the blocks of the multiband scattering matrix S under an 
uncoupled hole regime, acquiring a more simple form: 



ck*a — f3*f3 
(3*S-ja* 
(38* ^YoL 
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(6) 



where N' = N/2. Consequently, transmission and reflec- 
tion amplitudes have reduced expressions respect to 

t = l/cx* with an = a/f.. + , , 

r = -5^^j with (5(7)ji = (5(7)fl„ + z(5(7)/^. 

Rewriting the phase time ([S]) for this specific case, we ob- 
tain: 

h , da J.. do.Rii 
a/. . ■ 



dE 
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(8) 



Expressions like ([7]) were previously reported for the reso- 
nant tunneling and band mixing in multichannel superlat- 
tices of propagating electron systems [ll| . 

3. Numerical results and discussion: tunneling time 
and giant conductance 

For the numerical calculations that will be presented 
here we shall consider our layered media with barrier thick- 
ness Lb = (10 - 180) A, and well width L,„ = (50, 150) A. 
For the barrier height we will have either 14 — 0.498 eV or 
Vb ~ 0.550 eV . The cases where the in-plane wavenumber 
(parallel to the interfaces) ~ 0, with vanishing band 
mixing effects, and the cases where 7^ 0, with coupling 
phenomena effects on the quantum transport (3ll44. 45 1 will 
be separately study. Focusing to study the behavior of the 
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Figure 2: Phase time r^i for a tunneling of holes through a single 
barrier of = 20 A, as a function of the hole incident energy at 
Kt = with 14 = 0.498 eV. The solid (dashed) line stands for 
i = fe/i-)-3/2 ("i— 1/2) direct path. In the broken interval a negative 
dip extents up to 0.4 ps. 
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Figure 3: Direct phase time th for hole tunneling through a single 
barrier as a function of the hole incident energy for several values of 
the barrier thickness at Kt ~ 0, with Vb = 0.498 eV and i = 1. 



phase time tu and the free motion time = nlcm'^/fi{k,)h 
(being Ic the cell longitude, m,* the effective mass of holes 
and {kz)h the eigenvalues of the scattering problem solved 
with the MSA[3[) as a function of the incoming energy 
Ei, the barriers height and thickness were fixed at Vb = 
0.23 eV and Lf, = 30 A. respectively. Both, a double bar- 
rier {n = 2 cells) and a superlattice of ri = 8 cells were 
considered. For consistency with the reported results for 
electrons a 30% concentration of Al was fixed inside 
the barrier: {GaAs/Alo.sGao.rAs/GaAs)"^. 

3 A. Uncoupled hole regime phenomena 

In Figure [21 we plot the phase time for Ih and hh 
through a single potential barrier. It is clear from this 
figure that the Ih modes (dashed line) are faster than hh 
modes (solid line), although the difference is not large. In 
an effort to validate our numerical simulation with experi- 
mental data, we compare the results of Figure [5] with those 
of Dragoman et. al. [3| The tunneling times measured for 
holes, for null or small magnetic field in a ballistic-regimen 
nano-device, are of the order of 10"^'^ s. To reproduce 
that situation, we canceled the interaction between holes 
by assuming Kt ~ (uncoupled regimen), in the absence 
of magnetic field. In our calculation the phase time for 
hh and Ih is also of the order of lO^^'^ s, which is in 
an acceptable qualitative agreement with experiment and 
it is just one order less than the traveling time for elec- 
trons through a single barrier [4^, as was detected exper- 
imentally for a DBRT [i^. We have observed broaden 



peaks of the phase time for /1/1+3/2 direct transition at 
E K, 0.7 eV and for lh_i/2 direct transition at E Ki 0.9 eV. 
We have observed, although not plotted, that the trans- 
mission coefficient maxima are located at comparable en- 
ergies. It was shown in Ref.(ll| that, the tunneling time 
follows the band structure. This feature of hole tunneling 
agrees also with calculations of Ref.j^^l, where a particu- 
lar case of the envelope function approximation was used. 
Evaluating the tunneling times for energies in the range 
0.40 eV !^ E !^ 0.60 eV we have found, ahhough not 
shown, a reduction of tu in 3 orders of magnitude with re- 
spect to that shown in Figure[5]when the in-plane momen- 
tum changes as Ak-t- ^ 0.001 A. This behavior of r.,, was 
earlier detected in experiments for a hole passage through 
a slightly different system [4^ -a double asymmetric quan- 
tum well (QW)-, and is considered to arise from valence 
subbands mixing (48| . 

Figure [3] displays a raising of the phase time while Lh 
increases in the case of uncoupled regime {kt ~ 0) of holes 
incidence. The same is achieved for light holes, although 
not shown in this figure. This behavior agrees qualitatively 
with the results of an experiment on time-resolved lumi- 
nescence spectroscopy [43|. For E > 0.4 eV, tu is sensitive 
to the raising of and it is shown an oscillating behavior 
of the transmission phase time as function of the incident 
energy. The resonances shift to smaller energies, increasing 
their intensities with L^. It is assumed understood that, 
when the de Broglie wavelength of the hole propagating 
mode is a multiple of Lh, there is a resonance in the trans- 
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Figure 4: Phase time versus barrier thickness for the direct tran- 



sitions hh 



±3/2 



hh 



±3/2 



given as soUd diamonds/circles {E < 



Vt/E > VJ,), in the uncoupled case (kt ~ 0). A fragmented exponen- 
tial fitting is shown in dashed line (with full circles) a,t E = 0.551 eV. 
The barrier height Vt = 0.498 eV. 



mission and correspondingly the phase time delay raises 
at these energies. As Lf, increases, new resonances appear, 
thus the oscillating behavior is achieved. These oscillations 
correspond to the quasi-stationary levels of the virtual QW 
with increasing width. In that sense, from tu properties 
we are able to forecast a longer lifetime for the lower un- 
bound hole levels. On general grounds, we found this in 
correspondence with the idea of the quasi-stationary states 
lowering, whenever the virtual QW becomes wider. 

Figure 2] shows phase transmission time versus barrier 



thickness for the direct transitions hh 



±3/2 



hh 



±3/2 given 



as solid diamonds/circles {E <Vb/E > Vb), in the uncou- 
pled case. We prescribe appealing tunneling events of holes 
for E < Vb under uncoupled propagating regime. This fig- 
ure displays tu for resonant direct hhj.^/2 ~^ ^^±3/2 tran- 
sitions (diamonds), which rapidly becomes autonomous 
from barrier broadening as Hartman predicted for elec- 
trons [2^. These results are in accordance with the tran- 
sit time for optical pulses through opaque barriers (which 
transmit about 10~^ of the incident radiation) [13| and 
also with theoretical pre dictions for the delay time in one 
and two dimensions [58| . The circles in Figure H] corre- 
spond to numerical results within the MSA, for the tun- 
neling time of direct paths (indicated in the legend) for 
E > Vb- Here we must stress, that the resonant oscillat- 
ing behavior (on top of an average slope of smooth growth 
[59|), is an expected phenomena due to the barrier's inter- 
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Figure 5: Direct phase time tu for hole tunneling through a DBRT 
structure as a function of the energy of an impinging hole stream 
at Kt ^ 0, with 14 = 0.550, Lj, = 20 Aand Lu, = 50 A, for i = 
/i/i_l_3/2 (solid line), and i = lh_i/2 (dashed line). The inset displays 
transmission coefficients for the same direct transitions, through an 
identical DBRT heterostructure, as a function of the incident energy. 
The labelling for the hole levels follows that of the Reference |45| . 



ference effects with the continuum states, wide accepted 
as Ramsauer-Townsend oscillations For E > Vb we 

briefly investigated segments in Figure |31 pursuing a lo- 
cally exponential-like behavior of TjA, as was suggested in 
a reported experiment with holes [43|, within a rank com- 
parable with their measurements. The interval of barrier 
thickness length variance in their measurements, ranges 
from 30 A up to 80 A. Notice in Figure 2] a sectional fit- 
ting, using an exponential growth function (gH (dashed 
line with full circles). We do believe that oscillations are 
already intrinsically present in Heberle's, et al. accurate 
investigation [43| and could be shown by including more 
experimental points or enlarging in rank the search. 

Phase time through DBRT and superlattice 

Figure[5]is devoted to analyze briefly, the spectrum and 
the quasi-stationary bound states lifetime in a DBRT. The 
curves of Figure [S] (inset), show several transmission phase 
times Tu resonances for the uncoupled hole regime. In par- 
ticular we have analyzed a semiconductor heterostructure 
of the form GaAs/ {AlAs/GaAs)'^ /GaAs. It was shown in 
Ref. [llj that the tunneling time follows the band struc- 
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n = 2 



Table 1: Comparison to typical quasi-steady DBRT hole-level struc- 
ture with Vt = 0, 550 eV, = 10 A, and = 50 A. 

Levels [45] Resonances [eV] 

Experiment [53]. Theory (MSA). 



hhi 


0,028 


0,027 


Ihi 


0,073 


0,071 


hh2 


0,111 


0,110 


hhs 


0,245 


0,247 


lh2 


0,299 


0,307 


hh^ 


0,428 


0,434 



ture profile. Figure [5] (inset) displays clearly, &t E < Vt, 
the familiar discrete quasi-stationary hole-level spectrum 
of the embedded QW, widely accepted from hole quan- 



M Iii-i3- We have 



tum transport calculations [3|, \21 
found at low energies {E < Vb), that the th resonances are 
allocated at energies which correctly reproduce the hole 
spectrum previously reported for the DBRT experiment 
|53| , as expected (see Table [Ij . Notice the good quanti- 
tative agreement of our approach in Table [TJ with similar 
values reported elsewhere |53| . A previous compariso n [3|| 
with a robust and widely accepted theoretical model [27| . 
exhibits the advantages of the present approach. 

Tunneling time can not be quantitatively determined 



from the resonant tunneling line width, [57| nevertheless 
some qualitative prognosis is allowed. Notice the general 
trend of larger values for hh states phase time in compari- 
son to that of Ih states. This is in correspondence with the 
properties we have found for the transmission-coefficient 
resonances depicted in the inset of Figure [SJ which shows, 
as tendency, that hh resonances are thinner than the Ih 
ones. This observation foretells a stronger confinement for 
heavy holes, and we guess it could be used in resonant- 
tunneling diode design and nano-modelling j45| . 

Among the different formulations of tunneling time 
with recurrent experimental agreement |34[. the phase time 
have clearly advantages. In Reference [ll|, was demon- 
strated that the phase time associated with the passage 
(tunneling or not) of particles is a relevant quantity of 
the quantum transport, which was calculated and success- 
fully compared with precise experimental measurements 
[12 , 13 1 . Moreover in Reference [lH , some appealing effects 
for electron phase time were predicted, namely: superlu- 
minal phase propagation and resonant-band structure. In 
an attempt to distinguish similar behaviors for holes, we 
perform some numerical simulations, exhibited in the next 
two figures. 

In Figure [S] the phase transmission time T]^(2), through 
a single cell (n = 1) and a DBRT (n = 2) respectively, 
together with T2 are shown in connection with the free 
motion time Tf (dotted line). The last, is the propagation 
time of the quasi-particle in absence of scattering poten- 
tials. The direct paths under analysis, are indicated in 
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Figure 6: Top(Bottom) panel shows the evolution of the phase 
time, transmission probability, free motion time through n = 2 cells 
(DBRT) and the phase time through one single cell as function of the 
incoming energy for direct transitions of hh-j- 3/2 (^''^ 1/2)- With 
Vb = 0.23 eV, Lb = 30 Aand L„ 
Al was fixed inside the barrier. 



150 A. A 30% concentration of 



the legends of each panel. A resonant behavior for T2 is 
very clear for both hh states (top panel) and Ih states 
(bottom panel). In the forbidden energy region, T2 ap- 
proaches Ti (dashed red line), which represents its lower 
bound as had been foretold for electrons [ll|. The rela- 
tionship between T2 and T2 (blue line) resonances - as one 
can see from ([H])-, justify T2 to accurately reproduce, the 
quasi-stationary states of the embedded GaAs QW. An- 
other appealing feature arises for forbidden energy regions 
{T2 vanishes), where T2 < Tf. The opposite is predicted 
for allowed energies. For hh states (top panel) of incom- 
ing energies around 100 meV, (r/ — T2) ~ 0.4 ps , while 
for Ih states (bottom panel) (t/ — T2) ~ 0.1 ps. This ear- 
lier arrival time for hh and Ih, was predicted for electron 
tunneling through a same heterostructure (ll| . and sug- 
gest a more speedily passage of holes through a DBRT 
like {GaAs/Alo,3Gao,rAs/GaAs)'^. 

It remains to notice that {Tf)j^j^ > {Tf)^^ for the DBRT, 
having a bigger difference at low energies, with respect to 
that of high energies. Another remarkable distinction is 
the width of T2 resonances, which are noticeable widest 
for Ih. This difference yields to a larger confinement for 
hh inside the embedded QW, with respect to that for Ih. 
Thus greater lifetime for hh in the scattering region is 
expected. 

In Figure [7] the phase time Tg and transmission coef- 
ficient Tgj through a super lattice (n = 8) undoubtedly 



6 



A 2 



± 




1^^1/2 "'^1/2 



A 1 




0.2 0.3 

E [eV] 



Figure 7: Top (Bottom) panel shows the evolution of the phase 
time, transmission probability, free motion time through n = 8 cells 
(DBRT) and the phase time through one single cell as function of the 



incoming energy for direct transitions of hhj^ 3^2 ("iq: 



Vt = 0.23 eV, Li, = 30 Aand L„ 
Al was fixed inside the barrier. 



1/2J 



With 



150 A. A 30% concentration of 



follow a resonant-mini-band structure. Direct hole transi- 
tions were considered and appear indicated in the legends 
of each panel. For finite periodic systems -which is the 
case envisioned here-, a finite number of intra- band en- 
ergy levels [ill are typically observed, at variance with 
the continuous band structures predicted for infinite peri- 
odic systems. In the mini-gaps, Tg vanishes and rg tends 
to Ti, which is similar to the DBRT case discussed before. 
On the other hand, in the mini-bands region rg exhibits 
a resonant behavior approaching r/, which behaves as a 
lower bound [lH. In the scattering regions, rg increases as 
the incoming energy approximates to any allowed energy 
level. While, in the forbidden regions the holes propaga- 
tion is away faster, reflecting the evanescent behavior and 
the lack of hospitality of the barrier. These evidences for 
hh and Ih phase time, resemble those for electrons as was 
foretold for the tunneling time in a similar superlattice 

ini. 

Figure [S] exhibits two dissimilar behaviors of th for di- 
rect transitions hh±^/2 ~^ ^^±3/2 as function of the barrier 
thickness. For energies under the potential barrier (full tri- 
angles) Tii becomes independent of the barrier thickness, 
being consistent with Hartman's classical premonition for 
electrons [2^. In the DBRT, th becomes autonomous 
starting from if, « 60 A, meanwhile, in a superlattice 
of 71 = 4 cells, this quantity saturates at about L}, « 50 A. 
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Figure 8: Phase transmission time as a function of the barrier thick- 
ness for heavy hole transition ^''±3/2 within the uncou- 
pled regime, kt ~ 0. Both incident energies E: less (full triangles) 
and greater (empty triangles) than the barrier height Vt = 0.498 eV 
were fixed at 0.475 eV and 0.551 eV, respectively. The two bottom 
and top panels represents the numerical values of this time for the 
DBRT and the superlattice of n = 4 cells respectively. 



We can also note a decrease of the phase time in the region 
(10 - 15) A for the DBRT and around (10 - 17) A for the 
superlattice, showing the attractive effect of the barrier, 
as a reminiscence of the similar behavior widely quoted in 
the literature for electrons through a simple barrier jl, 
For energies higher than the potential barrier (empty tri- 
angles ) Tii, in the DBRT and in the superlattice, exhibits 
a resonant oscillating behavior similar to those of the sin- 
gle barrier (see Figure U]) . This is a consequence of an 
interference effect between the reflected wave and the con- 
tinuum states, known as Ramsauer-Townsend oscillations 
as commented above. By comparing Figure 3] and Fig- 
ure m it is clear that oscillations [1, [60| remain robust for 
the DBRT. The robustness of this feature vanishes for the 
superlattice, provides that well-resolved oscillating curves 
are no longer observed for n = 4 cells. 
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Giant conductance phenomena 



-output channel hh^^j2 
■ output channel //i_i/2 




kt = 1 * 10-'' A- 
Vh = 0.498 eV 



2.0 



E (eV) 



Figure 9: The bottom panel plots the 4 outgoing channel conduc- 
tances Gj through a single barrier heterostructure with Lj, = 20 A, 
where i = 1,4 (solid line) represents the hh±^/2 channels, respec- 
tively, and j = 2, 3 (dashed line) represents the channels, 
respectively. The top panel displays the cosine of phase shift for the 
same transitions. 

It is quite well established in the literature the strong 
relation between the transmission magnitudes and the wave 
function's phase [13, 0, E^l- Specifically, the giant con- 
ductance events are nearly related with phase shifts and 
phase interference phenomena [l^. It has been shown 
for various mesoscopic systems, the correlation between 
the variations of the conductance and an even function 
of the phase shift, due to an alternatively appearance of 
maximum-minimum valued sharp peaks or dips, whenever 
an even functional of the phase shift, changes by almost 
TT radians 



17, \m 54, 55, m 



Here, the cosine of the 
phase shift was chosen to verify, whether strike maximum- 
transmission peaks occur or not, in a negligible conduc- 
tance background. Commonly, the cos{Oij) is plotted, as 
it is a suitable even function to show evidences of giant 
conductance events [13, [l3 • 

Figure [§] (bottom panel) , presents the conductance Gj 
to direct channels — > hh^;^/2 (solid line), //i±i/2 
lh±i/2 (dashed line), and the cosine of phase shift for the 
same transitions (top panel) for the uncoupled hole regime 
as function of the incident energy. It is clear that the max- 
imum values in the conductance for both output channels 
are always accompanied by an extreme change in the co- 



sine of the phase difference. One can see, that the curve 
for hh output channel (solid line) changes two times in tt 
radians and in the bottom panel appear two Gj maximum 
correspondingly (at E « 0.72 eV and at £^ « 1.25 eV). 
Meanwhile for the Ih output channel (dashed line), there is 
only one change (at E sa 0.98 eV). Direct paths contribu- 
tions are the largest ones, as it follows from vanishing band 
mixing effects w 0. The maximum values shown 

in the bottom panel of Figure ^ for an_ even function of 
A6 (cosA6'), at even multiples of tt 



i7,ha, M 54,^ 



are not enough to confirm an abrupt increment of trans- 
mission. Indeed, there is nothing particularly large about 
hole transmission through a single cell, as that reported 
for electrons in finite periodic systems with time reversal 
symmetry [l^. Besides, in the single-cell case, does not 
exist a resonant mechanism that allows hole conductance 
to become giant through the barrier. However, it is tempt- 
ing to have a fairly reference frame to compare to what we 
will see next for a DBRT. Worthwhile noticing -although 
not shown here-, the largest contribution of the imaginary 
scattering amplitude component, that give rise Gj to max- 
imize in the same energy intervals, as we have shown in 
the Figure |9] bottom panel. 



K = 0.498 eV 
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Figure 10: Dependency of the 4 outgoing channels conductance and 
the cosine of phase shift with hole incident energy, through a DBRT 
and a superlattice with n = 4 cells of width Lj, = 20 A, and height 
Vt = 0.498 eV. 

Well-established time reversal symmetry of the (4 x 4) 
Kohn-Liittinger model, is not fulfilled by the (2 x 2) sub- 
spaces as was recently demonstrated [56|. Events of GG 
have been observed in many finite periodic systems, which 
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possess time reversal invariance but as yet not observed in 
time asymmetric ones. Probably, this explains why other 
accurate calculations made in the (2 x 2) Kohn-Luttinger 
subspaces, did not report GG events for holes [43, 4^ 47, 
[Hot . Figure [TO] displays the conductance and phase shift 
behavior as a function of the incident energy, for direct 
transitions to hh±^/2 (solid line) and lh±i/2 (dashed line) 
output channels. The Gi curves exhibit evidences pointing 
out the existence of well-defined resonances for E < Vb. It 
is relevant in this case, the presence of striking conduc- 
tance events, as can be clearly noticed in the panels for 
n = 2 of Figure \W[ Indeed, peaks lh2 and /1/14, have 
found to have a clear signature of GG phenomenon since: 
(i) they are sharp resonances, (ii) they keep a simple rela- 
tion with an even function cos{0ij)n=2 of the transmission 
phase {Oij)n=2, thereby peaks' locations at energy, are co- 
incident with those of {0ij)n=2 whenever the last changes 
by almost tt radians in a relative short energy interval [ITI - 
[19I . 23 . [6^ . and finally (iii) there is a strong phase in- 



terference mechanism between propagatin g or evanescent 
modes, giving rise to resonant transmission [l9j. As follows 
from the scattering theory, the GG mechanism in a DBRT, 
derives from the interference between states, whenever the 
incoming modes energy perfectly matches with one of the 
quasi-stationary states within the embedded QW of the 
scatterer system. This agreement produces the essential 
conditions for energy to maximize, i.e. a resonant trans- 
mission occurs. Sharp peaks become visible at energies 
E ^ 0.29 eV (dashed line) for direct transitions of Ih, and 
for E = 0.42 eV (solid line) for direct transitions of hh. 
We have disregarded the Ihi sharp resonances and oth- 
ers in the low energies region, because they do not reach 
maximum and also the corresponding phase shifts differ 
significantly from tt radians. In addition, it is necessary 
to mention the loss of the GG phenomenon in the super- 
lattice. As one can see in the top panels of Figure [TO] (for 
n = 4 and E < Vb), the {9ij)n=i changes are not about tt 
radians. In this case the leak of maximum transmission for 
the superlattice becomes visible. We conjecture the exis- 
tence of a competitive mechanism when the number of cells 
grow, and consequently a phase coherence effect rises. For 
higher energies, the direct Gi might give rise to a plateau, 
and we foretell the probable existence of another related 
striking feature: the giant resistance. Despite a more am- 
biguous interference process in this sharp dips are 
expected as a consequence of a transmission rate decrease, 
due to a resonant transfer of holes flux toward evanescent 
modes. 

It is likely to underline that the energies of the reso- 
nances -whether they are giant or not-, accurately repro- 
duce the quasi-bounded hole spectrum as was published 
for the DBRT [23|, were another approaching frame was 
applied. This is to say that, none GG phenomenologi- 
cal analysis can possibly be complete if a resonant means 
is missing. In this sense, to endorse the GG evidences 
showed in Figure [TO] the most immediate tool within the 
MSA, is the TM technique to estimate the hole states en- 



ergies at the QW of the DBRT. The basic idea here, is 
to search several states of the embedded finite quantum 
well (fQW), by cautious tuning those of the external in- 
finite quantum well (iQW). The last we had achieved, by 
imposing that Zl,r slowly tends to zi_4, respectively (see 
Figure [TT|) . This method approaches the iQw boundary 
conditions to those of the fQW, and is a reasonably good 
sketch for an estimation of some fQW levels. 
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Figure 11: Schematic representation of the quasi-bound-state levels 
(horizontal lines) and valence band lineup, within an infinite walls 
quantum well. The ground state wave function (dotted line) is repre- 
sented qualitatively. The energy is taken as positive for convenience. 

Figure [TT] shows schematically, holes quasi-bound-state 
levels (horizontal lines) and valence band lineup, within 
the iQW. For a numerical estimation of the quasi-stationary 
levels within the fQW, we consider a cell enclosed between 
the zl and zqs points. The barriers and well have the same 
dimensions of the Figure [TO] Thus, it is undemanding and 
straightforward to reproduce a DBRT or a superlattice, by 
a periodic concatenations of cells from zoa [3||. The latter 
considers facts of equivalence for points zi{z2) and 23(^4), 
respectively, as they have the same materials to the left 
and to the right. As boundary condition for the iQW, is 
standard to impose the wave function goes to zero at the 
points (zl ~ Az) and (z^ -|- Az) 



(y9(zL - Az) = (^(zr + Az) 



0, 



(9) 



being Az the splitting distance between these points and 
the infinite potential. In the framework of the TM ap- 
proach, the latter requirement can be written as: 





M/d(^;R + Az,ZL- Az) 





V'(z) 



(10) 



(zl-Az) 



Therefore, the non trivial solutions of the expression (|TU| 
can be numerically determined, by solving the equation 



Det{MKu{zn + Az, Zl - Az)} = 0, 



(11) 



being TWad the amplitude- derivative matrix block of the 
TM Mfd- Up to this end, it is turn clear from Table [2] 
a reasonably good agreement between our estimation for 
the energies of several quasi-steady levels of the fQW and 
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Figure 12: Two-probe conductance through the four direct accessible 
channels as a function of the incoming energy on the kt ~ limit, 
for a superlattices of n = 7, 8, 9 cells (from bottom to top) in the 
range of energies between [0.2 — 0.45] eV. 



the resonant-band structure displayed in Figure 1101 Now 
it is undoubted, the very existence of these fQW levels, 
and therefore their role as a fair resonant mechanism of 
the GG phenomena in the DBRT we claim for. The lh2 
and hhi energy values showed here, are a truthful prove 
of the GG phenomena, predicted for the uncoupled hole 
transport in the DBRT for sharp resonances allocated at 
similar energies (see Figure [TU]). 

Figure [T^] presents the dependency of the conductance 
G of hh and Ih through a superlattices of ri = 7, 8, 9 cells 
with the energy of the incident flux. The holes travel 
along all four accessible direct channels hh±^j2 /i/i±3/2 



Table 2: Quasi-steady fQW hole-level for a DBRT structure, quoted 
by Eq.^TTJ, with Vi = 0.498 eV , Lt = 20 A, and = 50 A. 



Levels (Labeling from Ref.[45]) 


Resonances [eV] 


hh2 


0.1362 


lh2 


0,2978 


hhi 


0,4112 



a 

H 

a 
o 
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Figure 13: Phase transmission time through one single cell versus 
barrier thickness for the direct and crossed passages of holes in a 
coupled case (kt = 5 X 10"* A"!) for E = 0.475 eV < Vb = 
0.498 eV. 



(solid line) and Z/i±i/2 — >■ lh±i/2 



12 (dotted-dashed line). It is 
wide open visible the thinner peaks of the hh resonances, 
respect to those of the Ih ones, since we are restricted 
on the uncoupled regime and thereby ft,ft,-spectrum states 
are the closest due to their larger effective mass. Fur- 
ther outstanding behavior, is the filter-like effect that can 
be observed on both flavors of holes in the selected re- 
gion. Indeed, the barriers become opaque for the hh/lh 
states at low/high energies, leading to a stronger confine- 
ment for the hh/lh at low/high energies. This fact yields a 
larger/shorter lifetime for hh/lh inside the allowed region 
of the scatterer system for low energy, at variance with the 
opposite for higher energies. Dealing with a superlattice of 
n = 9 cells (topmost panel), we found an hh — Ih overlap 
region at the vicinity of 0.3 eV. The phenomenology de- 
scribed in Figure I12[ allows us hypothesize it as concerns 
applications in hole-based electronics, pursuing maximum- 
minimum conductivity response as function of the charge 
carrier effective mass. 

3.2. Coupled hole regime phenomena 

Let us considers first the Hartman effect under this 
regime. Figure[T3]depicts the phase transmission time as a 
function of the barrier thickness for several crossed transi- 
tions paths, namely: hh^^/2 ~^ ^^1+1/2: '^-1/2 ~^ ^'1+1/2: 



lh_ 



1/2 



hh. 



3/2 



and hh 



+3/2 



lh_ 



1/2: 



of the sixteen 



possible ones, through a single-cell scatterer system. This 
plot shows the phase transmission time curves as a func- 
tion of the barrier thickness Lb, given fixed E <Vb. Notice 
the changes of Tij through each path while Lb increases, 
until Tij reaches saturation values. This autonomous be- 
havior of the phase time with Lb, unambiguously gives 
rise to the paradoxical Hartman effect [29| , which remains 
robust, even for a relative strong coupling of holes. Ow- 
ing to simplicity we dropped the remain transitions, pro- 
vided they follow a similar tendency. Worthwhile un- 
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Figure 14: The bottom panel shows the crossed-path contributions to 
the conductance through the first outgoing channel, Gi , for a DBRT 
with barriers of width = 20 A, and height Vt = 0.498 eV. The 
top panel displays the cosine of two phase shifts that match with the 




same outgoing channel, hh^^^2 ~ 
both in the same range of energy. 
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Figure 15: The bottom panel shows the the crossed-path contri- 
butions to the conductance through the second outgoing channel, 
G2, for a DBRT with barriers of width L;, = 20 A, and height 
Vj, = 0.498 eV. The top panel displays the cosine of two phase 



shifts that match with the same outgoing channel, lh_ 



and lh_i/2 hh^ 



1/2 



■ Ih. 



3/2 



both in the same range of energy. 



1/2 



derlying some characteristics of transient interval for 
several paths. Crossed transitions hh+z/2 — >■ lh+1/2 (+) 



and lh_ 



1/2 



(up empty triangles), show negative 



We turn now to discuss how robust are the hole gi- 
ant conductance features, in the presence of valence-band 
mixing. Figure \TM and Figure [T51 present the conductance 



slope before reach saturation. Meanwhile the crossed path and the cos(%)„=2 through the outgoing channels /1/1+3/2 
^^-1-3/2 ^ ^'^-1/2 (down full triangles), exhibits a double 
change of slope. These behaviors clearly departs (in the 
transient interval) from the classical Hartman premoni- 
tion for electrons and deserve more investigation to shed 
light on this topic. In the presence of interference effects 
(kt 7^ 0) at w 0.55 eV, we have found abrupt tunneling 
phase changes when A9 w tt radians, which seems to be 
an anomalous behavior of the phase time. In Figure [T^ 
we have found an overlap of Tij for several crossed tran- 
sitions. For example: with j : 1(2) — > 3 : i at Lf, « 20 
A, also with j : 1(2) 2(3) : i at Lb « [110, 140] A, and 
with J : 1 — 7> 2(3) : i at Li, ~ [70,80] A. In some cases 
the effective mass changes, and the incoming channels be- 
long to a so called Kramer- degenerated hole states 
while the outgoing channels are Kramer- /oiu ones. It can 
be readily observed, similar Tij matches, for other crossed 
transitions at different values of Lf,. These coincidences 
relate to the strong hh — Ih propagating modes mixing, to 
the simultaneous treatment of the 4-coupled channels of 
the present approach, and to hh ^ Ih interaction with an 
effective potential due the barrier broadening. Ac- 
cordingly, whenever a match in takes place, the time 
spent by the hh and Ih in the classically forbidden region, 
shall be the same. 



and lh_i/2, respectively, as a function of the incoming en- 
ergy under interference between channels. We had taken 
Kt = 10"'^ A"-'^. Despite the presence of sharp peaks at 
E < Vb (qualitatively similar to those showed in Figure 
[TU)) . as one can observe at E ^ 0.16 eV, E = 0.22 eV 
and E = 0.24 eV in FigfHl it is not seen a neat concor- 
dance with the widely accepted GG formulation posted 
above [13, [H, [H, [H, Ell , and illustrated in FigUHl Firstly 
there is not a nice correspondence between the highest 
peak of Gi at the vicinity of £; = 0.24 eV (see FiglTil 
bottom panel) and the fluctuation of (%)„=2, which for 
this peak is sensitive less than tt radians (see FiglH] top 
panel). On the other hand, the values of Gi,2 exhibited 
in Fig[T3] and FigfTSl respectively, are in average almost 
one order smaller than those shown by FigflUl being the 
comparative behavior observed in Fig llSl even more restric- 
tive in this sense, with respect to FiglTU] Following these 
evidences, we conjecture so far, that the GG events are 
not straightforward signatures for hole tunneling travers- 
ing a DBRT, under a selected range of parameters in the 
coupled-particle regime. 

Additional remarkable differences become clear from a 
comparison of this trend in the coupled hole regime and 
Figure [TOl For instance, there is a small shift between 
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Figure 16: Full contributions to the conductance through the second 
outgoing channel, G2 , for a DBRT with barriers of width Lf, = 20 
A, and height Vt = 0.498 eV. 



the position in energies of the peaks and the changes in 



the cos( 



^2- We guess that the last reflects the loss 



of maximum transmission when the valence-band mixing, 
i.e. the interference between paths is incremented. Be- 
sides, the Gn oscillations spectra in the low-energy inter- 
val shown in Fig ll5l seem to be of the Fano-type leading 
to assume the existence of interference of embedded-QW 
discrete levels with a continuum. The later gives rise to 
characteristically asymmetric peaks in transmission spec- 
tra [H. 

The porpoise of Figure [T71 is to briefly analyze fur- 
ther anomalous features in the transmission phase time 
spectra. Fig 1 171 displavs in panel (a) the cos{9ij)n=i for 
the crossed Ih — hh transition — > ^^-3/2 (dashed 
line, i : 3 — 4 : j). In panel (a) with dotted line, we 
plot the difference of the the same magnitude for the path 
i : 3 — >■ 3 : j with respect to the path i : 3 -> 4 : j. 
The bottom panel (b) plots tunneling phase time for the 
crossed-transition i : 3 — >■ 4 : j through a single barrier 
heterostructure. In the presence of interference effects 
(kt ^ Q) aX E k, 0.55 eV, the phase changes abruptly 
when « vr radians, as can be seen from the dashed 
line in the top panel of Figure [T71 Both curve in panel (a) 
show a clear correspondence of abrupt changes in cos 9ij 
with observed peaks(dips) of T ij, which seem to be an 
anomalous behavior of the phase time, as they resemble 
Fano-like resonances due theirs asymmetry hallmark (63 |. 
In the neighborhood of 0.8 eV [see panel (b)] we have found 
Fano-like positive (negative) peaks (dips) of with i j, 
displayed in the bottom panel of Figure [TTj They are as- 
sociated with changes of the phase in almost 7r/2(7r/3) 
radians, and are plotted in the top panel with dashed line. 
We also have observed, although will not be shown here, 
local increments of transmission coefficients in the energy 
vicinity, where negative dips appeared, though regretfully 
they need more investigation to derive proper conclusions. 

Most interestingly, though, we found indications that 
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Figure 17: The bottom panel plots tunneling phase time for the 
crossed-transition j : 3 — >■ 4 : j through a single barrier heterostruc- 
ture at Kt = 0.001 with L;, = 20 A and Vf, = 0.498 eV. 
Top panel shows the cosine of phase shift for the same transition 
J : 3 — > 4 : i (dashed line) where the first number stands for the 
input channel j and the other stands for the output one i. In dotted 
line the same magnitude is shown for the direct path j : 3 — > 3 : i 
with respect to the crossed path j : 3 — > 4 : i. 



might be connected to possible giant conductance events in 
the single-cell case, which we omit graphically here because 
of brevity. There were found two prompt shifts in trans- 
mission phase amplitude, equivalent to an even multiple of 
TT (|A0.y | = 67r) at « 0.01 A~^, similar to those plot- 



ted in Figure [TUl One of them, appear in Ih^ 



1/2 



lh_ 



1/2 



intra-band transition (z : 3 — >■ 2 : j) when Li, varies within 
[10 — 20] A range. The second one, was observed through 
^^-3/2 ~^ ^^+1/2 inter-band transition (i : 1 — > 4 : j) when 
Lb varies within [120 — 130] A range. Then, maxima of Gj 
could be expected for those crossed paths, not only vary- 
ing the incident energy, but the barrier thickness as well. 
However, a clear resonant mechanism in this case remains 
a puzzle and further investigations need to be addressed to 
clarify wether the effect of giant conductance for holes, is 
expected to arise under valence-band coupling phenomena 
during transport trough a single-cell system. 
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4. Concluding Remarks 

Relevant and striking quantum transport features of 
holes, directly connected to the phase, the phase shift 
and the phase resonant interference phenomena are: the 
giant conductance-resistance events, the Hartman effect, 
and the earlier phase propagation in the barrier. Sev- 
eral results nicely resemble those previously reported for 
electrons while some predictions of the MSA model 
are in acceptable qualitative compliance with the experi- 
ments HI, m, llll and with theoretical predictions [3 62|. 
Theoretical evidences of the giant conductance phenom- 
ena for decoupled hole transmission through DBRT were 
exhibited. In this case a trustworthy resonant mechanism 
supporting this anomaly, is the alignment of the incident 
propagation energy with one of the quasi-bond hole states 
of the embedded quantum well in the DBRT. This mech- 
anism was confirmed with an independent numerical cal- 
culations of several quasi-bound levels. The giant con- 
ductance effect vanishes in the presence of valence-band 
particles coupling and by tuning the number of superlat- 
tice layers as well. The hole phase transmission time for 
a {GaAs / Alo,3Gao,7As / GaAs)"' superlattice, shows a res- 
onant mini-bands structure in the scattering regions, hav- 
ing the free motion as an inferior limit time. Moreover, in 
the forbidden regions the holes propagation became away 
faster and was shown tu to be inferiorly bounded by the 
single-cell phase time ri. In the DBRT mini-gaps regions, 
Tii moves forward from the free motion time, suggesting a 
more speedily phase propagation within the barrier. The 
Hartman's classic prediction is observed robust for null 
and finite hole band mixing. Was presented an alternate- 
selective confinement strength independently for both fla- 
vors of holes, implying filter-like effects on the effective 
mass via the manipulation of the incident energy in the 
uncoupled hole regime. 
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